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^ ■ Abstract 

\Q ■ We describe some interesting effects observed during the evolution of nonequi- 

librium systems, using domain growth and glassy systems as examples. We 

y—t \ breafly discuss the analytical tools that have been recently used to study the 

dynamics of these systems. We mainly concentrate on one of the results ob- 
■ tained from this study, the violation of the fluctuation-dissipation theorem 

CN ■ and we discuss, in particular, its relation to the definition and measurement 

^ ! of effective temperatures out of equilibrium. 

ON : 

One of the major challenges in physics is to understand the behaviour of systems that 
are far from equilibrium. These systems are ubiquitous in nature. Some examples are 
phase separation, systems undergoing domain growth, all types of glasses, turbulent flows, 
! systems driven by non-potential forces, etc. All these systems are "large" in the sense that 
they are composed of many, N —> 00, dynamic degrees of freedom. Apart from succeeding in 
predicting the time evolution of their macroscopic properties, one would like to know which, 
if any, of the thermodynamic notions apply to these nonequilibrium cases. 

Systems undergoing domain growth, or phase separation, provide the best known exam- 
ple of a nonequilibrium evolution.0 Take for instance a magnetic system with ferromagnetic 
interactions in contact with a thermal bath. If the bath temperature is very high the sample 
is in its paramagnetic phase and the magnetic moments, or spins, point in random directions. 
If one next cools down the bath, and hence the sample, through a transition temperature 
T c , the system enters the low temperature phase and starts forming domains or islands of 
the two ordered phases, say up and down. For defmiteness, let us fix the final temperature 
to be < T < T c . At any time t w after crossing the transition at the initial time, two 
types of dynamics appear: (i) fast fluctuations of some spins, due to thermal fluctuations, 
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FIG. 1. Three two-dimensional slices of a domain-growth system at three different waiting 
times t Q < t w i < t w2 < tw3- See the text for a discussion. 



inside the otherwise fully ordered domains; (ii) slow motion of the domain walls leading to 
the growth of the averaged domain size L(t w ). If the size of the sample is infinite, in real 
life very large, the nonequilibrium domain-growth process can take so long that the sample 
simply does not equilibrate in the time-window that is accessible experimentally. In other 
words, below the critical temperature T c one always has t obs < r EQ with r 0BS the observation 
time and r EQ the equilibration time. The two types of dynamics itemized above are clearly 
seen in Fig. 1 where three two-dimensional slices of a system undergoing domain-growth 
are displayed. The pictures are obtained at increasing waiting times after the quench. One 
sees the domains growing as well as the existence, in each of the snapshots, of some reversed 
spins inside the otherwise ordered domains. 

Scaling arguments have been extensively used to describe the dynamics below T c ; they 
are based on the assumption (sometimes derivation) of the evolution of the averaged domain 
size L(t w ) and on further proposals for the space and time-dependence of the correlation 
functions. 

Understanding the physics of glassy materials is perhaps a problem of intermediate 
difficulty.!'! Glassy materials can be of very many different types; one has for instance 
structural glasses,! orientational glasses,! spin-glasses,!!! plastics,! gels and clays,! glycerol]! 
etc. Their hallmark is that below some transition range they also fall out of equilibrium. 

The easiest way of preparing a glassy system is again through an annealing.! This is 
implemented by decreasing the temperature of the bath with a given cooling-rate. Take 
for example the case of a molecular liquid. For high enough bath temperatures the sample 
is in its liquid phase and it achieves equilibrium with the bath. At an intermediate bath- 
temperature range the liquid avoids the crystallization transition, enters a metastable phase 
and becomes a super-cooled liquid, that is to say a liquid with some peculiar properties as, for 
example, a extremely high viscousity. At an even lower bath-temperature the liquid cannot 
follow the pace of the annealing and falls out of equilibrium, it becomes a glass. If one stops 
the annealing at any temperature below this range the system stays in its glassy phase for 
practical purposes forever and is typically an amorphous solid. We talk about a "transition 
range" since the transition might not be clearcut but depend on the cooling-rate. Actually 
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one can form a glass of probably any substance by choosing a fast enough cooling-rate. 
Many other routes to the glassy phase are also possible. 

There have been proposals to describe the evolution of some glasses, notably spin-glasses, 
with scaling arguments based on domain growth ideas.00 The assumption is that the glassy 
dynamics is simply given by the growth of domains of two competing ground states. However, 
it has been very difficult to prove (or disprove!) either experimentally or numerically that 
this is indeed the scenario: no "ordered structures" have been identified in general as the 
growing phases. 

Thus, domain growth, phase separation and glassy materials are all "self-sustained"[] out 
of equilibrium systems. If one follows their time-evolution, keeping all parameters fixed, in 
particular the bath-temperature, some of the main features observed during their nonequi- 
librium evolution are: 

Slow dynamics. The evolution is very slow. "One-time quantities", as the energy-density, 
approach their asymptotic limit with some slow decaying function, say power law, logarith- 
mic or more complicated. It is very important to notice though that even if these one-time 
quantities can get very close to their asymptotic values, this does not mean that the systems 
get frozen in a metastable state: they are not equilibrated in a restricted region of phase 
space characterised by these asympotic values. This is most clearly demonstrated by the 
measurement of "two-time quantities". 

Two-time quantities and physical aging. The measurement of these quantities prove that, 
even if one-time quantities approach a limit, the system is still changing in an important 
way. 

One can distinguish two types of two-time quantities. Those measured during the free 
evolution of the system, quantifying the spontaneous fluctuations, such as any two-time 
correlation function, and those measured after applying a small perturbation, such as any 
response function. These quantities depend on both times involved in the measurement and 
not only on the time-difference. This shows that the systems neither are equilibrated with 
the bath nor have approached equilibrium in any metastable state. They are indeed rather 
far from equilibrium. 

The measurement of the spontaneous fluctuations is quite easy to implement in a nu- 
merical simulation.lll One prepares the sample at an initial time t a and lets it evolve until a 
waiting time t w when the system configuration is recorded. One then lets the sample further 
evolve and computes, at all subsequent times t = r + t w , the correlation function between 
the reference configuration at t w and the configurations at r + t w . These curves depend on 
both t w and r and they are not invariant under time translations showing that the system is 
out of equilibrium. Furthermore, the decay as a function of r is slower the longer t w . This 
is the phenomenon called physical aging: the younger (older) the sample the faster (slower) 
the decay. 

The result of the measurement of a local auto-correlation function is very easy to vi- 
sualize for a domain growth. Take again the case of ferromagnetic domain growth, the 
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dynamic variables are the Ising spins that we encode in a time-dependent iV-dimensional 
vector 4>{t) = (4>i(t), . . . , 4>N(t)) (the index i = 1,...,N labels the spins) and the local 
auto-correlation function is just the scalar product of two configurations evaluated a dif- 
ferent times, NC(t,t w ) = (<fi(t) ■ 4>{t w )). For Ising spins, the auto-correlation function is 
normalized to one at equal times. A departure from one measures how different are any 
two configurations as those shown in Fig. 1. For any fixed t w the auto-correlation has two 
distinct regimes depending on the time-difference r = t — t w . Let us choose a waiting time 
t wl and plot C(r + t wl , t wl ) as a function of r. The curve has a first fast decay from one to a 
bath-temperature dependent value qEA(T) (the Edwards- Anderson temperature-dependent 
parameter). This corresponds to the decorrelation associated to the fast flipping of the 
spins inside the domains. In this regime r is small compared to an increasing function of 
the domain size g(L(t w )). When r increases and becomes of the order of g(L(t w )) one starts 
seeing the motion of the domain-walls, i.e. the growth of the domains, and the decay slows 
down. If one repeats this calculation choosing a longer waiting-time t w2 > t wl , and its 
associated reference configuration, one observes that the first decay is identical to the one 
for t wl though it lasts for longer, and that the second regime is notably slower that the one 
associated to t wl . These features can be easily understood. While r is smaller than g(L(t w )) 
the dynamics takes place only inside the domains as thermal fluctuations. The domain walls 
are ignored and the correlation behaves as if the system were a patchwork of the two equi- 
librium states. The correlation function decay is then independent of the waiting-time and 
approaches qEA{T) = m^ q {T). However, after a time-difference of the order of g(L(t w )) the 
system realizes it has domain walls, the subsequent decay is associated to the motion of 
the walls and is nonequilibrium in nature. The decay gets slower the longer the t w simply 
because the size of the domains reached at t w is larger. 

For glassy systems the correlation functions have exactly the same qualitative behaviour 
though, as already mentioned, it is not easy to decide if there is any type of order growing. 
Plots like the one displayed in Fig. 2 have been obtained for an impressive number of 
glassy models of different nature. Some of them are the 3D Edwards-Anderson modeljl! 
a polymer meltjl! a polymer in a random potential,0 a binary Lennard- Jones mixturejl! 
etc. Furthemore, the kind of curves were found in several sandpile models and other kind of 
systems.lll 

The measurement of "dc-response" functions or, more precisely, integrated dc-responses 
is what is usually done experimentally. The starting procedure is similar to the prece- 
dent: one prepares the sample at an initial time and lets it freely evolve until t w . At this 
waiting time one applies a small, constant, perturbation and then measures the associated 
integrated response of the system as a function of r for different t w . For example, exper- 
imenting with spin-glasses one applies a small magnetic field and measures the increase of 
magnetization,! manipulating with polymer glasses one applies a stress and measure the 
tensile creep compliance,^ etc. In all cases, the integrated responses are studied as functions 
of t w and t and they all show aging effects that manifest in a similar way as in the corre- 
lation measurement. There is a first increase of the time- integrated susceptibility towards 
a value Xi(T) that does not depend on the waiting-time while there is a second increase of 
the time-integrated susceptibility towards the equilibrium value Xeq(^) that is waiting-time 
dependent. 

This can again be simply visualize in the domain growth problem. The first regime 
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FIG. 2. The local auto-correlation function for wait-times t w \< t w2 < t w3 < t w4 < t w5 < t w $. 
See the text for a discussion of the different regimes. 



corresponds to the response of the spins inside the domains, i.e. to the response of the 
full system taken to be roughly a patchwork of independent equilibrium states. The second 
regime instead corresponds to the response of the domain walls. Since their density decreases 
as time elapses, one expects this nonequilibrium response to vanish in the long witing-time 
limit. 

"Ac-response" measurements are also usually performed experiment allyffl In these ex- 
periments one applies an ac small field of fixed frequency u at the initial time to and keeps 
it applied until the measuring time t w . The relation with the previous results is given by 
identifying u ~ 1/r. The out of equilibrium character of the evolution is given by an explicit 
t w dependence in the relaxation of the in and out of phase susceptibilities. 

It is important to notice that these effects are physical aging as opposed to chemical 
aging. Physical aging is totally reversible: it suffices to heat the sample above the transition 
range and cool it back again below it to recover a fully rejuvenated system. 

The comparison of spontaneous fluctuation, e.g. a correlation function, to induced fluctu- 
ations, measured by its associated response function, is well established for systems evolving 
in equilibrium. Indeed, this relation involves the temperature of the bath and it is called 
the fluctuation- dissipation theorem. However, for systems that are far from equilibrium this 
relation does not necessarily hold. 

The problem of reaching a theoretical understanding of nonequilibrium physics is im- 
portant both from a practical and a theoretical point of view. It is obvious that for some 
applications, one would like to predict the time-evolution of the samples with great precision 
and avoid undesired changes that depend on the, sometimes unknown, age of the samples. 
From an analytic point of view, domain-growth and glassy materials are a pattern of out 
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of equilibrium systems whose properties one could try to capture with simple models or 
simplified approaches to more complex models. The predictions thus obtained can then be 
experimentally (or numerically) tested in real systems. Importantly enough, one could also 
try to extend some of these predictions to other nonequilibrium systems such as those ex- 
ternally driven. Some connections between glassy systems and driven systems are discussed 
in Ref. 0. 

How can one modelise domain growth or glassy systems? The "microscopic" constituents 
and interactions and, consequently, the microscopic models differ from glass to glass. Spin- 
glasses are composed of magnetic impurities (spins) that occupy fixed random positions in 
the sample and interact via RKKY interactions; polymer glasses are composed of polymers 
(strings) that interact via potentials of (oversimplifying) Lennard- Jones type. In the former 
case, there is quenched (time- independent) disorder in the system that is associated to the 
random positions of the spins, that give rise to random interactions between them (the 
RKKY interactions oscillate very rapidly with the distance between the magnetic impurities 
and change sign in an almost random manner). In the latter case nothing can be interpreted 
as being quenched disorder. Though one could expect that these two systems (and other type 
of glasses) behave very differently, the experimental results as well as the recent theoretical 
developments show that their dynamical behaviour is indeed rather similar. In other words, 
one can identify certain quantities that have the same qualitative behaviour. 

In order to describe the dynamic evolution of a classical system in contact with an 
environment one starts by identifying the relevant variables of the system that evolve in time. 
One then proposes a Langevin equation, with noise and friction mimicking the coupling of the 
system to the thermal bath, to determine the time-evolution of the time-dependent variables. 
These variables are the (soft) spins in the case of a spin-glass, the monomer positions in the 
case of the polymer glass, etc. This procedure leads to a set of N, the number of dynamic 
variables in the system, coupled differential equations of second (if inertia is included) or 
first (if inertia is neglected) order. Obviously, this huge system cannot be solved and one 
has to resort to some alternative method to further advance in the analysis. 

Indeed, one would like to obtain information about "macroscopic" quantities, that can 
be related to experimental measurements, instead of following the erratic motion of any 
microcopic variable. The two-point correlation or response functions are macroscopic quan- 
tities we would like to monitor. A well-known theoretical method, known under the name 
of Martin-Siggia-Rose (MSR) formalismjll allows us to obtain a generating functional, and 
from it the Schwinger- Dyson integro-differential equations, for these quantities. This method 
can be applied in complete generality, to any classical model. For models with non-linear in- 
teractions of finite range it recquires the calculation of an infinite series of diagrams that one 
cannot in general resum and express in terms of two-point functions in an explicit formJi^l 
One then faces the problem of choosing some approximation scheme to simplify this series 
expansion. 

Before discussing how to deal with this problem, let us describe the formalism used to 
study a quantum system in contact with an environment. The Schwinger-Keldysh closed- 
time path (CTP) formalism was developed to monitor the nonequilibrium time-evolution of 
a quantum system, and to obtain information about two-time quantities.^ The environment 
is usually modelized by a set of harmonic oscillators (infinitely many for each variable in 
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the system) with a spectral distribution of frequencies.^! The coupling of system and bath 
is usually chosen to be linear but of course more general situations can be considered. In 
this way, one obtains the CTP generating functional that, as in the classical case, involves 
a series expansion that, in general, cannot be obtained explicitly. (The classical limit of the 
CTP generating functional is the MSR generating functional.) 

Typically, two routes are followed to approximate these generating functionals. They are 
equally applicable in the classical and quantum case and are the following:! 

• The microscopic models, namely the starting Hamiltonians, are simplified in such a 
way that the construction can be carried through and that explicit equations can be 
derived. This is the choice made when one uses, for example, the large N limit of a 
0(iV)-model to describe domain growth, fully connected spin models to describe spin- 
glasses or when one embeds a finite dimensional manifold in an infinite dimensional 
space to describe an interface motion, the motion of a polymer in a random medium. 

• The microscopic models are realistic but some approximation scheme is chosen to 
select, from the infinite series, a still infinite subset of diagrams that can be resummed 
to yield an explicit set of dynamic equations for correlations and responses. Many such 
recipes exist in the literature,!!^ some of them are the mode-coupling approximation, 
the direct interaction approximation, the self-consistent screening approximation, etc. 

These two procedures yield the same "form" of coupled integro-differential causal equa- 
tions. Actually, in some cases one can show that a simplified microscopic model with infinite 
range interactions (e.g. the p spin-glass model) yields the same dynamic equations that an 
approximation scheme (e.g. the mode coupling appoximation) applied to a more realistic 
model for a glassy material.il@ The structure of these equations is always the same: there 
might be a second time-derivative term if there is inertia, some terms describe the interac- 
tion of the system with the bath and some other integral terms describe the interactions in 
the system (through the self-energy and vertex). It is the explicit form of the self-energy 
and vertex that is selected by the model or the approximation. 

Once one has the equations governing the evolution of the two-time quantities, for all 
values of the parameters, the question then arises as to which is the phenomenology that 
they describe. 

A combination of analytic and numeric methods are used to study these equations. One 
can attempt a numerical solution taking advantage of the fact that they are causal. The 
solution shows that they present a dynamic phase transition at a temperature Td- Above 
Td, the solution reaches, in the long waiting-time limit, a stationary form. All two-time 
correlations and responses are functions of the time-difference only and are related through 
the fluctuation-dissipation theorem. The high-temperature dynamic equations for spin- 
glass were studied in detail by Sompolinsky and Zippelius0 for spin-glass models, by Gotze 
and collaborators!! for glass-models and the relation between these two was signalled and 
investigated by Kirkpatrick, Thirumalai and Wolynesi! in a series of beautiful papers. 

Below Td, a drawback of the numerical method is that, due to the slowness of the 
dynamics and the memory of the system, one cannot reach very long time intervals. The 
numerical solution gives us hints about the structure of the solution but does not give us 
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extremely precise information about more detailed features such as the two-time scaling laws, 
etc. Nevertheless, the numerical solution sufficed to show that below two-time functions 
start depending on the waiting-time and that aging is captured by these equations.il 

Below T<i, and in the asymptotic limit of long waiting-time, an analytical solution was 
developed first for the p spin-glass mo dell and later for other mean-field disordered models 
such as Sherrington-KirkpatrickEl or the motion of manifolds in infinite dimensional random 
potentials.^!! 

One of the main ingredients of this solutionis concerns the fluctuation- dissipation the- 
orem that relates, in equilibrium, the spontaneous to the induced fluctuations. Indeed, if 
one follows the dynamics of a classical system that is in equilibrium with a bath, one can 
easily show that 
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±JL{0{t)0{1?)) 0{t-H) = ~C{t,H) 9(t-t>), (1) 



n=0 



with 0(t) any observable taken to have zero mean for simplicity and h an infinitesimal field 
acting a time t' that modifies the energy of the system according to V — > V — hO and that 
is not correlated with the equilibrium configuration of the system. 

In the glassy phase, this relation does not hold. This does not come as a surprise since 
the equilibrium condition under which it can be proven does not apply. What really comes 
as a surprise is that the modification of the relation between response and correlation takes 
a rather simple form for domain-growth and glassy systems. 

A way to quantify the modification of FDT in the out of equilibrium phase and to use 
it to classify different systems is the following. Let us integrate the response function over 
a time-interval going from a waiting-time t w to a final time t: 

X {t,t w ) = ( alt'R{t,t') . (2) 



This yields a time-integrated susceptibility that is exactly what is measured experimen- 
tally. Next, we compare this integrated-susceptibility to the auto-correlation function. In 
equilibrium, one can use FDT to show that 

X (t,t w ) = ^(C(t,t)-C(t,t w )) . (3) 

Hence, if one draws a plot of x against C, for increasing t w , using r = t — t w as a pa- 
rameter, in the large t w limit the plot will approach a straight line of slope — 1/T join- 
ing (lim^oo C(t, t), 0) and (0, % EQ ). From now on and without loss of generality we take 
limt^oo C(t,t) = 1. Any departure from this straight line signals a modification of FDT and 
a departure from equilibrium. 

The analytic solution of simplifed models shows that, in the nonequilibrium phase, this 
construction converges to a limiting curve given by 

t w —*oo,c(t,t w )=c ^ ^ 7 1 EFF (C) 
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where T EFF (C) is a function of the correlation C . We shall discuss the notation and justify 
the name of this function below. In the large t w limit two distinct regimes develop in the x 
vs C curve. There is a first straight line of slope — 1/T, joining (1,0) and (^^(T), xi(T)). 
This characterises what is called the FDT regime. The straight line then breaks and the x 
vs C curve goes on in a different manner. The subsequent behaviour depends on the model. 
Indeed, three families have been identified: 

• Models describing domain growth like, for example, the O(N) model in D dimensions 
in the large N limit. In this case, one follows the local correlation NC(t, t w ) = ((Mx, £)• 
<p(x, t w )) and its associated local susceptibility. The plot for C < qEA(T) is flat." The 
susceptibility gets stuck at its value XiCO while the correlation continues decreasing 
towards zero. The same result holds for the Ohta-Jasnow-Kawasaki approximation to 
the A0 4 model of phase separation.!^ 

• Models describing structural glasses like, for example, the so-called F p _i models of the 
mode-coupling approach or the p spin-glass models. In this case the x vs C plot, for 
C < Qea{T), is a straight line of slope larger than — 1/T.H 

• Models describing spin-glasses like, for example, the Sherrington-Kirkpatrick model. 
In this case the x vs C plot, for C < qEA(T), is a non-trivial curve.0 

This "classification" in three families has been checked numerically for more realistic 
models. Many numerical simulations using either Montecarlo (MC) techniques or molecular 
dynamics (MD) have shown that several models fall into the expected cathegories. Some 
models belongingto the first group are the 2D Ising model with conserved and non-conserved 
order parameterjl the site diluted ferromagnet and the random field Ising modelH and the 
2D Ising model with ferromagnetic exchange and antiferromagnetic dipolar interactions.il 
The binary Lennard- Jones mixtures are a standard model for theglass transtion. Both MC 
and MD simulations show that they belong to the second class." Finally, MC simulations 
of finite dimensional spin-glass models, the three and four dimensional Edwards- Anderson 
model, yield the third kind of behaviour. Another particularly interesting problem, relevant 
for the physics of dirty superconductors,^! is the one of a manifold diffusing in a random 
potential. The analytic prediction using an infinite dimensional embedding space depends 
on the nature of the quenched random potential, namely on it being short or long range 
correlated.^ This prediction is partially confirmed by the simulations in finite dimensional 
transverse space with the proviso of a very interesting modification that is not captured 
by the infinite dimensional approach.0 Besides, numerical simulations of lattice-gas modes 
with kinetic constraints^ and sandpile models^! also show FDT violations. 

Once this modification of FDT in the nonequilibrium situation is identified, several ques- 
tions arise, all connected with the initial purpose of checking which thermodynamic concepts 
can be applied, perhaps after some modifications, to the nonequilibrium case. In the follow- 
ing we discuss three interesting issues. 

• Why is there always a two-time regime, when C first decays from its equal times value 
to qEA(T), where FDT holds? 
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For the domain-growth problem the presence of this piece is easy to justify. In this 
time-scale one only sees the dynamics and the effect of the perturbation inside the domains. 
Since one can then ignore the presence of domain walls, the equilibrium relation between 
correlation and response is expected to hold. Of course, one cannot easily extend this 
argument to a more general situation. There is however a totally general reason for having 
FDT when C > qea{T) and it is the following. 

For any system in contact with an environment, with bounded correlation functions and 
without non-potential forces^ the departure from FDT is bounded by@ 

\T X (t,t w )-C(t,t) + C(t,t w )\<Kj t yt' (--L^l) 1 (5) 

where K is a finite constant and 7 the friction coefficient that characterises the coupling to 
the bath. The Kubo 7i-function is defined as@ 

H(t') = Jdc/>d^P(ct>,^t)(T\nP( ( j>^,t) + V(<j>) + ^^ , (6) 

with P(<p, 0, t) the time-dependent probability distribution, V(<p) the potential energy and 
m a mass. The 7i-function satifies Ti. < for all times and it vanishes only for the canonical 
distribution. 

From this bound one sees that if 7~t(t) falls to zero faster than 1/t no FDT violations are 
allowed in the long t w limit since the right-hand-side in Eq. (|5]) vanishes. Instead, if TC(t) 
falls to zero in a slower manner, FDT is imposed by the bound for small time-differences 
but violations are allowed for longer time-differences. This argument proves that there is 
always a region of correlations close to C = 1 in which FDT holds, even for a system that 
is not close to equilibrium. 

• Can one identify the slope of the plot with an inverse effective temperature and call it 
-1/T EFF (MJ = -l/T EFF (C)? 

About ten years ago, in the context of weak-turbulence, Hohenberg and Shraiman0 
proposed to define an effective temperature through the departure from FDT. However, a 
detailed analysis of this quantity and its properties was not given in this reference. 

Indeed, one expects that any quantity to be defined as a nonequilibrium effective tem- 
perature must fulfill the requirements associated to the intuitive idea of temperature. The 
first property to check is if this effective temperature is measurable by a thermometer that 
is weakly coupled to the system, in a statistical manner, at any chosen waiting time." This 
property can be proven by studying the time-evolution of the thermometer coupled to M 
identical copies of the system, all of age t w , and by verifying that this equation becomes a 
Langevin equation in the presence of a thermal bath characterised by a coloured noise with 
correlation given by the system's correlation and response given by the system's response. 



2 Other bounds can be found if diffusion and/or non-potential forces are allowed.@El 
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Thus, if the system has several time-scales characterized by different values of the effective 
temperatures^ 

C(t, t w ) = C FDT (t, t w ) + (t, t w ) + (t, t w ) + . . . (7) 
R(t,t w ) = R FDT (t,t w ) + R^(t,t w ) + R^(t,t w ) + . . . (8) 

with 

R FDT (t,t w ) = ±JL C Fm (t,t w )9(t-t w ) R {i \t,t w ) = -LJLc®(t,t w )9(t-t w ) (9) 
T ot w Tw ot w 

one can select which value is measured by choosing the internal time-scale of the ther- 
mometer. Say, for example, that the thermometer is a harmonic oscillator of internal fre- 
quency uj . Then, one chooses the system time-scale to be explored, and hence the value of 
the effective temperature to be measured, by comparing u to t w . 

Many desirable "thermodynamic" properties of T EFF defined in this way can also be 
checked,c§ for example: 

(i) T EFF controls the direction of heat flows. 

(ii) T EFF controls partial equilibrations between observables in a system that evolve in the 
same time-scales and interact strongly enough. 

(iii) Let us take two different glasses, in contact with a single bath of temperature T. These 
glasses are constructed in such a way that when they are not in contact each of them 
has a piecewise T EFF (C) of the form 

rpSYSTl/^-YN _ / T if C > q" E \ rp S YST2(r<\ — i T if C > (7£A /lfV\ 

FFF (C) " \ T« if C < q%\ eff (C) - j T(2) {iC< ^ (10) 

with T( 2 \ One can then reproduce the experiment of setting two observables 

in contact by coupling these two systems through a small linear coupling between 
their microscopic variables. The result is that above a critical (though small) value of 
the coupling strength the two values the effective temperatures below qEA equal while 
below the same critical value of the coupling strength the values remain unaltered. 
One concludes that if the two observables interact strongly the systems arrange their 
time-scales in such a way to partially thermalise. 

The presence of non trivial effective temperatures in glycerol out of equilibrium is 
presently being checked experimentally by Grigera and IsraeloffS Their results show 
that, at fixed measuring frequency u ~ 8 Hz, this system has an effective temperature 
T EFF > T = 180K until measuring times of at least 10 5 sec, that is to say of the order of 
days! (Note that the bath temperature T is below T c = 187 K.) 



See Ref. [30] for a precise definition of two-time scales. 
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Further support to the notion of effective temperatures comes from the study of the effect 
of quantum fluctuations on the same family of models. Below a critical line, that separates 
glassy from equilibrium phases, and in the slow dynamic regime, one finds violations of the 
quantum fluctuation dissipation theorem. These are characterised by the replacement of 
the bath temperature by an effective temperature T EFF (t,t w ). The effective temperature is 
again piecewise. It coincides with the bath-temperature T when C is larger than qEA and 
it is different when C goes below qEA- This nonequilibrium value has the nice property of 
being non-zero even at zero bath-temperature. Again, this result can be interpreted within 
the domain growth example. Whenever one looks at short time- differences with respect to 
the waiting-time one explores the quantum and thermal fluctuation in the bulk, i.e. one 
observes a quantum equilibrium dynamics that satisfies the quantum FDT. Instead, when 
t is comparable to g(L(t w )) one observes the domain wall motion. These are macroscopic 
objects for which quantum fluctuation do not have a strong effect. This can be seen, for 
example, in the form of the FDT violations: they look classical though with an effective 
temperature that depends on the strength of quantum fluctuations. 

• Do effective temperatures in out of equilibrium systems emerge from a symmetry 
breaking? 

In the classical case, one can study the structure of time-scales and effective temperatures 
with the help of the supersymmetric formulation of stochastic pro cesses. 0H Indeed, it is 
well-known that the effective action in the MSR generating functional is invariant under a 
supersymmetric group (with a possible symmetry breaking due to the initial condition). 

In the kind of glassy systems we deal with, there is a neat separation of time-scales in 
the long waiting-time limit. This allows us to separate the dynamics in the fast scale from 
the dynamics in the slow time-scales. The equation governing the slow time-scales, have an 
enlarged symmetry: they acquire an invariance under super-reparametrizations. The only 
solution that respects the large symmetry is a trivial, constant one. Hence, in order to 
have non-trivial dynamics in the long waiting-time limit, the system has to spontaneously 
break the super-reparametrization invariance. One can prove that the choice of effective 
temperatures is intimately related to the spontaneous breaking of this invariance.il 

A similar analysis in the quantum case remains to be developed. 

In conclusion, we have summarized some interesting features of the slow out of equilib- 
rium dynamics of domain growth and glassy systems. We have explained why these features 
arise in the domain-growth case. A similar understanding has not been reached for glassy 
systems yet. With the purpose of developing a "visual" understanding of glassy physics, 
a careful analysis of the statistics and organisation of the configurations visited by a glass 
model during its nonequilibrium evolution is in order. 

The use of simplified models or, alternatively, self-consistent approximations to more 
realistic ones have yielded a number of very interesting results and new predicitons. In 
particular, these models capture much of the aging phenomenology of glassy systems. Sur- 
prisingly enough, even puzzling effects of temperature cyclings during aging in spin-glasses, 
and the absence of these effects in other kind of glasses, can be described by fully-connected 
modelsB Some of these new predictions, notably the modification of FDT, have been tested 
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numerically and experiments are now being performed. Obviously, it is desirable to go 
beyond these approximations and study more realistic models in finite dimensions. This, 
however, is a very difficult task. 

There have been innumerable attempts to define a temperature for an out of equilib- 
rium system. In particular, in the context of glassy materials, a "Active temperature" is 
often introduced to describe some of the experimental findings.^ The effective temperature 
discussed in this article has the most welcome property of being measurable, hence being 
open to experimental tests. As far as we have checked the definition, it also has the welcome 
property of conforming to the common prejudices one has of a temperature. Of course there 
are still many open questions related to it. Just to mention one, let us say that it would 
be very interesting to extend the analytical experiment of "coupling a thermometer to a 
system" to the quantum case. 
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